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The atomic structure of icosahedral B4C boron carbide is determined by comparing existing 
infra-red absorption and Raman diffusion measurements with the predictions of accurate ab initio 
lattice-dynamical calculations performed for different structural models, a task presently beyond 
X-ray and neutron diffraction abihty. By examining the inter- and intra-icosahedral contributions 
to the stiffness we show that, contrary to recent conjectures, intra-icosahedral bonds are harder. 



Covalently bonded solids based on boron, carbon, or 
nitrogen form the hardest materials presently known [Q , 
and B4C comes third after diamond and cubic BN, with 
the advantages of being easily synthesized and stable up 
to very high temperatures . Hence it is used as abrasive 
or shielding material sustaining extreme conditions, while 
^"^B-enhanced ceramics are used in nuclear reactors. 
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FIG. f . Atomic structure of B4C. 



The atomic structure responsible for these properties 
is rather unique an arrangement of distorted BuC 
icosahedra located at the nodes of a rhombohedral Bra- 
vais lattice [Rim space group, Fig. |l|). This is a mod- 
ification of the a-Bi2 structure, which accommodates 
three-atom chains along the (111) rhombohedral axis, 
linking different BuC icosahedra. Two inequivalent crys- 
tallographic sites exist in the icosahedron. The six atoms 
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which form the top and bottom triangular faces of the 
icosahedron sit at the polar sites, and are directly linked 
to atoms in neighboring icosahedra; the other six corners 
of the icosahedron form a puckered hexagon in the plane 
perpendicular to the [111] axis, and their symmetry- 
equivalent sites are called equatorial. Each one of the six 
equatorial atoms is linked to an inter-icosahedral chain. 

An experimental determination of the atomic struc- 
ture of B4C is still lacking |Q,pf . Neutron diffraction [HJ^] 
cannot distinguish ^"'^B from C because their scattering 
lengths are too close Q . X-ray diffraction has allowed to 
identify a C-B-C chain |^,^, but the location of the re- 
maining C atom in the BuC icosahedron remains unset- 
tled because the X-ray form factors of boron and carbon 
atoms are also too close. A recent attempt to resolve the 
atomic structure of B4C using photo-emission and X-ray 
absorption spectra |^ did not yield conclusive results. 

In this paper the determination of the atomic struc- 
ture of B4C is addressed by a combination of first- 
principles calculations, based on density-functional the- 
ory (DFT) and density-functional perturbation theory 
(DFPT), with infra-red absorption and Raman diffusion 
results. We have considered three ordered configurations 
of B4C: in the chain configuration, all the icosahedra are 
entirely constituted by boron atoms — B12 — while all the 
carbon atoms lie entirely in the inter-icosahedral chains — 
C-C-C; in the polar configuration, the chains are as- 
sumed to be C-B-C, in agreement with X-ray diffraction 
data , while one of the polar atoms of the icosahe- 
dra is substituted by a carbon atom — BuC; the equa- 
torial configuration is similar to the polar one, but in 
this case the substitution involves an equatorial atom. 
These configurations are highly idealized, and some de- 
gree of substitutional disorder is expected to occur in ac- 
tual samples. Nevertheless, their vibrational properties 
are sufficiently different so as to allow one to discriminate 
between different structural models by comparing the ob- 
served Raman and infra-red spectra with the predictions 
of accurate lattice-dynamical calculations. We finally de- 
termine the equation of state of B4C and point out that 
the recently postulated inverted molecular behavior ||ic| ] 
is not consistent with our theoretical findings. 

Calculations were performed within DFT and DFPT 
Iprf , using the local-density approximation and the plane- 
wave pseudo-potential method. The pseudo-potentials of 
boron and carbon were respectively the same as in Ref. 
1^ and (isj . Plane- waves up to a kinetic energy cutoff of 
65 Ry have been included in the basis set. The irreducible 
wedge of the Brillouin zone has been sampled with 10 and 
2 points for static and dynamic properties respectively in 
the chain model, and 20 and 3 points in the polar and 
equatorial ones. The structural parameters reported in 
table I have been obtained by minimizing the crystal en- 
ergy with respect to the size, shape and internal degrees 
of freedom of the unit cell. The chain model has -R3m 
symmetry, while the substitution of one B atom in the 
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TABLE I. Comparison between the structural parameters 
of B4C as calculated for three structural models [Chain, Po- 
lar, Equatorial: see text) with those observed experimentally. 
a is the lattice parameter (A), a is the angle between the 
rhombohedral lattice vectors (degrees), while &„ indicates the 
length of the n-th bond (A) (see Fig. Ill) . 
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icosahedron induces a small monoclinic distortion, which, 
for the polar and equatorial configurations, amounts to 
1.8% and 0.5% of the rhombohedral cell length respec- 
tively, and to 1% and 0.1% of the rhombohedral angle 
respectively. The calculated data agree well with X-ray 
1^ and neutron diffraction Q data (see Table I) for all of 
the three configurations. An agreement better than 3% 
is found for almost all of the bond lengths, which depend 
very little on the configuration. Exceptions to this are 
the intra-chain bond (bond #1 in Fig. [^) whose length 
would be underestimated by 10% in the chain model and 
the chain-icosahedron bond (bond # 7 in Fig. |l]) which 
is predicted to be 5% too short both in the polar and the 
equatorial models. We interpret the failure to accurately 
predict the chain-icosahedron bond as due to our neglect 
of disorder effects in the chain. Occasionally, the chain 
could be B-C-C , a possibility which we did not con- 
sider. The calculated energies are also very close: those 
of B12-CCC and equatorial B4C are slightly higher by re- 
spectively 73 and 4 meV/atom as compared to polar B4C 
which would then result to be the ground state. However, 
our neglect of disorder and of finite-temperature effects 
does not allow us to draw any definite conclusions. 

In order to substantiate the hypothesis of a polar 
B4C, we have decided to compare the vibrational spec- 
trum predicted using state-of-the-art theoretical meth- 
ods (DFPT) with that observed experimentally. As the 
difference of the atomic masses of boron and carbon 
amounts to 11%, one expects to observe a difference on 
those vibrational modes which involve mainly one of the 
two sites. In panels a-c of Fig. |^ we display the in- 
frared spectrum as calculated for the three configurations 
and as deduced from reflectivity experiments |14| . The 
polarization-averaged absorption coefficient reads : 
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where Z*^ is the Born effective-charge tensor of s-th 
atom, c is the speed of hght, is the unit cell volume, a 
and P are Cartesian coordinates, j labels the vibrational 
mode, and loj , e^^ the corresponding eigen-frequency and 
displacement pattern, respectively. For the polar and 
equatorial configurations, the dynamical matrix has been 
averaged with respect to all the symmetry-equivalent 
substitution sites in the icosahedron, so that rhombo- 
hedral symmetry is recovered also in these cases, and 
{A2u) modes are infrared active when light is polarized 
perpendicular (parallel) to the [111] axis. Phonon life- 
time and other broadening effects have been empirically 
accounted for by replacing the 6 functions in Eq. (|^) with 
Lorentzians with a half width of 10 cm~^. The positions 
of the peaks agree only for the polar configuration. 
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FIG. 2. Panels a-c: infrared spectra of B4C as calculated 
for the three configurations considered in this work (solid 
lines) compared with that deduced from experiments p^ ] 
(dashed lines), (a) B12-CCC; (b) equatorial B4C; (c) polar 
B4C. Panel d: calculated spectrum of polar B4C fitted to ac- 
count for the experimental mixing of polarization. (dash) 
and A^u (dot-dash) contributions shifted for clarity. 



The agreement between the calculated and observed 
peak intensities is still rather poor, mainly due to our 
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TABLE II. Polar B4C: theoretical infrared active fre- 
quency Lo (cm~^) and relative integrated intensity Irei- The 
highest intensity of each polarization is normalized to 10. F 
is the half width of the Lorentzians used in Fig. 2d (cm~^). 
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lack of knowledge of the individual line-widths as well 
as of the polarization conditions of the experiment with 
respect to the orientation of the sample. A simple adjust- 
ment of the theoretical data turns out to bring them in 
good agreement with experiments. We have re-weighted 
all the peaks with respect to the ones by a same 
factor which is treated as a fitting parameter. Once this is 
done, we have kept the integrated intensity of each mode 
at its theoretical value, and we have utilized the widths 
of individual peaks as free parameters. The resulting 
values of the fitting parameters are reported in table II. 
The agreement between the experimental spectrum and 
the one predicted for the polar configuration is now very 
good (Fig. ||d), and only the peak at 950 cm~^ is missed 
by the theory. Its shape is asymmetric and its intensity 

1^ , hence suggesting 
that it is not a simple bulk lattice vibration. In fact, the- 
ory predicts an infrared mode at 972 cm^^ but its small 
intensity should hardly be observable (Tabic II) . It might 
be activated by surface effects and/or bulk defects. Ac- 
tually experiments are performed on hot-pressed micro- 
crystalline samples, where the defect density is very high. 
This leads to a degradation of the reflectivity spectrum 
with respect to the ideal (mono-crystalline) case, and the 
largest experimental error is expected at the LO minima. 

In order to further confirm the polar configuration as 
the dominant one in B4C, we now compare the Raman 
spectrum measured on a single crystal fl^ ] with the pre- 
dictions of our calculations (Fig. ||). Unfortunately, a 
complete quantitative comparison is difficult due to mode 
broadening, especially at high frequency. However, two 
modes at 481 and 534 cm^^ are rather sharp |l^, and 
they are properly accounted for only by the polar model, 
which predicts two peaks respectively at 498 and 543 
cm^^. The former is a rotation of the chain about an axis 
perpendicular to the [111] direction and the later is the 
librational mode of the icosahedron previously identified 
in a-boron | [T^ . These two modes have similar frequen- 
cies and widths, and they both involve angular atomic 
displacements, so that the errors made in the prediction 
of their positions are expected to be similar and to can- 
cel to a large extent. The distance of 54 cm^^ between 



varies from sample to sample |14 
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the two peaks is therefore a stringent test for the assign- 
ment, and only the polar configuration does account for 
it, while it is underestimated by a factor of two in the 
equatorial model (Fig. ^). At higher frequencies, struc- 
tural disorder is likely to be responsible for most of the 
observed broadening. Twinning or stacking faults — re- 
cently identified as pairs of neighboring twin defects — are 
observed in single crystals They distort the crys- 

tal only slightly because they preserve the rhombohedral 
symmetry and the structure of the icosahedron. In B4C 
ceramics obtained by hot pressing, however, these defects 
can be as close as a few lattice parameters , and the 
low-frequency part of the Raman spectrum is quite differ- 
ent from those of Fig. 3. In particular, while no Raman 
activity is expected in clean samples below 360 cm~^, 
a broad band systematically appears at low frequency in 
ceramics |17 1^, twice-peaked around 250 and 320 cm~^, 
with an intensity which differs from sample to sample of 
a same bulk composition |l7|. We attribute this band 
to a high density of states of disorder-activated acoustic 
phonons, in analogy to what is found, e.g., in ice VII [GOl. 
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FIG. 3. B4C Raman spectrum: (A) experiment (B) 
theory for (a) B12-CCC, (b) polar B4C, (c) equatorial B4C. 
Solid lines: Eg mode; dashed lines: Aig mode. 



Having shown that the vibrational spectra of B4C can 
be consistently interpreted in terms of a carbon atom oc- 
cupying one of the polar sites, we consider now its elastic 
properties, which have recently attracted some interest 
p0{ . Theoretical volume vs. pressure data - determined 
up to a an hydrostatic pressure of 10 GPa - have been 
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TABLE III. Theoretical and experimental bulk modulus 
(GPa) of the unit cell and of the icosahedron. Theoretical 
pressure derivatives are found to be: i3o=3.5. 
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fitted to a Murnaghan's equation of state: the resulting 
bulk modulus Bq is reported in Table III together with 
the 'compressibility of the icosahedron' evaluated as in 
Ref. . Our theoretical value is close both to ultrasonic 
measurements |2l|] and to unpublished X-ray data [ p2| . 
No evidence of the inverted molecular behavior reported 
in Ref. was found from our calculations. In Rcf. pO[ , 
a powder sample of NaCl and B4C was used. Since the 
two materials have widely different elastic constants, they 
might have experienced very different pressures if, as it 
is often the case, some strain is transmitted at the grain 
boundaries (Voigt model). As no correction for this ef- 
fect was made, the bulk modulus was probably found too 
low since at a given NaCl pressure (calibrant) the actual 
pressure experienced by B4C was higher. For this rea- 
son the apparent bulk modulus of the icosahedra, which 
was found smaller than Bq, might need to be verified in 
a quasi hydrostatic environment, with an accurate mea- 
surement of pressure. Moreover, twinning in boron allows 
the crystal to accommodate large strains and to prevent 
the formation of quasi-crystalline phases . We are in- 
clined to attribute the observed anti-molecular compres- 
sion to pressure inhomogeneities induced by the intensive 
twinning in B4C ceramics. An experiment on a single 
crystal would provide a definite answer to this issue. 

In conclusion, the atomic structure of B4C consists of 
one icosahedron BuC with the carbon atom staying at 
a polar site, and a chain C-B-C. B4C and a boron are 
not molecular nor — as we have demonstrated — inverted 
molecular crystals. Rather, they should be considered as 
members of a new class of covalently bonded materials. 

The authors are grateful to J. Loveday and D. Sime- 
one for the communication of unpublished results. N.V. 
thanks L. Zuppiroli for a fruitful discussion. 
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